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ABSTRACT 

We report the discovery of WTS-lb, the first extrasolar planet found by the WFCAM 
Transit Survey, which began observations at the 3.8-m United Kingdom Infrared Tele- 
scope. Light curves comprising almost 1200 epochs with a photometric precision of 
better than 1 per cent to J ~ 16 were constructed for ~ 60 000 stars and searched 
for periodic transit signals. For one of the most promising transiting candidates, high- 
resolution spectra taken at the Hobby-Eberly Telescope allowed us to estimate the 
spectroscopic parameters of the host star, a late-F main sequence dwarf (V=16.13) 
with possibly slightly subsolar metallicity, and to measure its radial velocity variations. 
The combined analysis of the light curves and spectroscopic data resulted in an orbital 
period of the substcllar companion of 3.35 days, a planetary mass of 4. 01±0. 35 Mj, 
and a planetary radius of 1.49l ; 18 Rj. WTS-lb has one of the largest radius anomalies 
among the known hot Jupiters in the mass range 3-5 Mj. 

Key words: Extrasolar planet, Hot- Jupiter, Radius anomaly 



* Based on observations collected at the 3.8-m United Kingdom Infrared Telescope (Hawaii, USA), the Hobby-Eberly Telescope 
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1 INTRODUCTION 

The existence of highly-irradiated, gas-giants planets orbit- 
ing within < 0.1 AU of their host stars, and the unexpected 
large radii of many of them, is an unresolve d problem in 
the t heory of planet formation and evolution (|Baraffe et al.l 
|2010T ). Their prominence amongst the 777 confirmed ex- 
oplanet^] is unsurprising; their large radii, large masses, 
and short orbital periods make them readily accessible to 
ground-based transit and radial velocity surveys, on account 
of the comparatively large flux variations and reflex motions 
that they cause. Exoplanet searches often focus on the detec- 
tion of small Earth-like exoplanets, but understanding the 
formation mechanism and evolution of the giant planets, 
particularly in the overlap mass regime with brown dwarfs, 
is a key question in astrophysics. 

This paper reports on WTS-lb, the first planet discov- 
ery from the WFCAM T ransit Survey (WTS; iKovacs et all 
120121 : iBirkbv et al.|[20lTh . The WTS is the only large-scale 
ground-based transit survey that operates at near infrared 
(NIR) wavelengths. The advantage of photometric monitor- 
ing at NIR wavelengths is an increased sensitivity to photons 
from M-dwarfs (M* < 0.6 Mq). These small stars undergo 
similar flux variations and reflex motions in the presence 
of an Earth-like companion as solar-type stars do with hot 
Jupiter companions. While a primary goal of the WTS is 
the detection of terrestrial planets around cool stars, WTS 
can also provide observational constraints on the mechanism 
for giant planet formation, by accurately measuring the hot 
Jupiter fraction for M-dwarfs. The light curve quality of the 
WTS is sufficie nt to reveal transiti ng super-Earths around 
mid-M dwarfs (|Kovacs et al.ll2012T ). but scaling-up means 
that any hot Jupiter companions to the ^80 000 FGK stars 
in the survey are also detectable. 

Theoretical models of isolated giant planets predict an 
almost constant radius for pure H+He objects in the mass 
range 0.5 — 10 Mj as a result of the equilibrium between 
the electron degeneracy in the core and the pressure sup- 
port in the external gas layers (Zapolsky fc Sarpeterlll969l ; 
lGuillodl20"05l ; ISeager et al.ll2007l ; baraffe et al.ll2010] ). Hence 
the larger radii of many hot Jupiters (hereafter HJ) must 
arise from other factors. Due to the proximity of these plan- 
ets to their host star, the irradiation of the surface of the 
planet is thought to play a major role in the so-called ra- 
dius anomaly, by altering the thermal equilibrium and de- 
laying the Kelvin-Helmholtz contrac tion of the planet from 
birth (e.g. IShowman fc GuiUotl 120021 ). This is supported by 
a correlation betwe en the mean planeta r y density and the 
incid ent stellar flux (lLaughlin et al ] |2010l ; lDemorv fc Seagerl 
l201ll ; lEnoch et alj|2012l ). However, it has been shown that 
this cannot be the on ly explanation for the radius anomaly 
(|Burrows et al.ll2007l ). and other sources must contribute to 
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the large amou nt of energy required to keep gas giants radii 
above ~1.2 Rj jBaraffe et alj|2010h . 

There are a number of physical mechanisms thought to 
be responsible for radius inflation, including (but not lim- 
ited to ): tidal heating due to the circularisation of close-i n 
orbits (IBodenheimer et al. I l200ll . 120031 ; IJackson et alll2008l ) , 
reduced heat loss due to enhanced opacities in the oute r 
layers of the planetary atmosphere (|Burrows et alj l2007h . 
double-di ffusive convection leading to slower heat trans- 
portation jChabrier et alj|2007l ; iLeconte fc Chabrierll2012l ). 
increased heating via Ohmic dissipation in which ionised 
atoms interact with the planetary magnetic field as they 
move along strong atmospheric winds (|Batygin fc Stevenson! 
I2010I ). and a slower cooling rate due to the mechanical 
greenhouse effect in which turbulent mixing drives a down- 
wards fl ux of heat simulating a m ore intense incident irra- 
diation |Youdin fc Mitchellll2010l ) . A radically d iffere nt ex- 
planation has been proposed by iMartfn et al.l (|201ll ) who 
point out a correlation between radius anomaly and tidal 
decay timescale and suggest that inflated HJs are actually 
young because they have recently formed as a result of bi- 
nary mergers. Studying the radius anomaly in higher mass 
HJs (> 3Rj) is useful as they are perhaps more resilient to 
atmospheric loss due to their larger Roche lobe radius. 

This paper is organised as follows: in Section 11.11 we 
briefly describe the strategy of the WFCAM Transit Sur- 
vey, while the photometric and spectroscopic observations 
are presented in Section 12.11 along with their related data 
reduction. Section [3] describes how we characterized the 
host star (Section 13. ip and determined the properties of the 
planet (Section 13. 2p . using a combination of low- and high- 
resolution spectra and photometric follow-up observations. 
A discussion on the nature and the peculiarity of this new 
planet, our conclusions and further considerations are given 
in Section^ Throughout this paper we refer to the planet as 
WTS-lb, and to the parent star and the whole star-planet 
system as WTS-1. 



1.1 The WFCAM Transit Survey 

The WTS is an on-going photometric monitoring cam- 
paign using the Wide Field Camera (WFCAM) on the 
United Kingdom Infrared Telescope (UKIRT) at Mauna 
Kea, Hawaii, and has been in operation since August 2007. 
UKIRT is a 3.8-m telescope, designed solely for NIR obser- 
vations and operated in queue-scheduled mode. The survey 
was awarded 200 nights of observing time, of which ~50% 
has been observed to-date, and runs primarily as a back-up 
program when observing conditions are not optimal for the 
main UKIDSS programs (seeing > 1 arcsec). The survey 
targets four 1.6 deg 2 fields, distributed seasonally in right 
ascension at 03, 07, 17 and 19 hours to allow year-round 
visibility. WTS-lb is located in the 19 hour field (hereafter 
19hr). The fields are close to the Galactic plane but have 
b > 5 degrees to avoid over-crowding and high reddening. 
The exact locations were chosen to maximise the number of 
M-dwarfs, while keeping giant contamination at a minimum 
and maintaining an E(B — V) < 0.1. Due to the back-up 
nature of the program, observations of a given field are ran- 
domly distributed throughout a given night, but on average 
occur in a one hour block at its beginning or end. Seasonal 
visibility also leaves long gaps when no observations are 
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Table 1. WFCAM J-band light curve. The full epoch list, which 
contains 1 182 entries, will be available electronically. 



Table 2. The host star WTS-1. 



HJD 


Normalized 


Error 


-2 400 000 


flux 




54317.8138166 


1.0034 


0.0057 


54317.8258565 


0.9971 


0.0056 


55896.7105702 


0.9979 


0.0064 



possible ijKovacs et al.ll2012r i. Since the WTS was primar- 
ily designed to find planets transiting M-dwarf stars, the 
observations are obtained in the J-band (~1.25/mi). This 
wavelength is near to the peak of the spectral energy dis- 
tribution (SED) of a typical M-dwarf. Hotter stars, with an 
emission peaking at shorter wavelength, are consequently 
fainter in this band. Interestingly, photometric monitoring 
at NIR wavelengths may have a further advantage as it is 
potentially less susceptible to th e effect of star-spot induced 
variability (|Goulding et al.ll20T2l ) . The discovery reported in 
this paper demonstrates that despite the survey optimisa- 
tion for M-dwarfs, its uneven epoch distribution and the 
increased difficulty in obtaining high-precision light curves 
from ground-based infrared detectors, it is still able to detect 
HJs around FGK stars. 



2 OBSERVATIONS 
2.1 Photometric data 

2.1.1 UKIRT/WFCAM J-band photometry 

The instrument used for this campaign is WFCAM 
|Casali et al.ll2007ft . which has a mosaic of four Rockwell 
Hawaii-II PACE infrared imaging 2048 x 2048 pixels detec- 
tors, covering 13.65' x 13.65' (0.4" /pixel) each. The detec- 
tors are placed in the four corners of a square (this pattern 
is called a paw-print) with a separation of 12.83' between 
the chips, corresponding to 94% of a chip width. Each of 
the four fields consists of eight pointings of the WFCAM 
paw-print, each one comprising a 9-point jitter pattern of 
10 second exposures each, and tiled to give uniform cover- 
age across the field. It takes 15 minutes to observe an entire 
WTS field (9 x 10s x 8 + overheads). In this way, the NIR 
light-curves have an average cadence of 4 data points per 
hour. 

A full description of our 2- D image processing and 
light curve generation is given bylKovacs et al.l ||2012t ) and 
closely follows that of llrwin et al. | (|2007T ). Briefly, a mod- 
ified version of the C ASU INT wide-field survey pipeline 
|lrwin fc Lewis! l200ll ) is used to remove the dark cur- 
rent and reset anomaly from the raw images, apply a flat- 
field correction using twilight flats, and to decurtain and 
sky subtract the final images. Astrometric and photometric 
calibr ation is performed us ing 2MASS stars in the field-of- 
view |Hodgkin et al.l |2009| ) . A master catalogue of source 
positions is generated from a stacked image of the 20 best 



Parameter 


Value 


RA a 


19h 35m 58.37s 


Dec a 


+36d 17m 25.17s 


l a 


+70.0140 deg 


b a 


+7.5486 deg 


fi a cos5 h 


-7.7+2.4 mas yr _1 


W b 


-2.8+2.4 mas yr" 1 



http: / /casu. ast.cam.ac.uk/surveys-projccts/wfcam/tcchnical/ 



a Epoch J2000; 

b Proper motion from SDSS. 



frames and co-located, variable aperture photometry is used 
to extract the light curves. We attempt to remove system- 
atic trends that may arise from fiat-fielding inaccuracies or 
varying differential atmospheric extinction across the image 
by fitting a 2-D quadratic polynomial to the flux residuals 
in each light curve as a function of the source's spatial po- 
sition on the chip. As a final correction, for each source we 
remove residual seeing-correlated effects by fitting a second- 
order polynomial to the correlation b etween its flux an d the 
stellar image FWHM per frame ( see llrwin et al.ll2007| for a 
discussion of the effectiveness of these techniques). 

The final J-band light curves have a photometric preci- 
sion of 1% down to J = 16 mag (~7% for J = 18 mag while 
the uncertainty per data point is just 3 mmag for the bright- 
est stars (saturation occurs between J ~ 12 — 13 mag). The 
unfolded J-band light curve for WTS-1 is shown in Figure [T] 
(all the measurements are reported in Table [TJ and its final 
out-of-transit RMS is 0.0064 (equivalent to 6.92 mmag). 



2.1.2 Transit detection algorithm 

We identified WTS-lb in the J-band light curves by us- 
ing the Box-Le ast-Squares trans i t sear ch algorithm OCCFIT, 
as described in lAigrain fc Irwinl (|2004h , which takes a maxi- 
mum likelihood approach to fitting generalized periodic step 
functions. Before inspecting transit candidates by eye, we 
employed several criteria to speed up the detection process. 
The first is a magnitude cut, in which we removed all sources 
fainter than J = 17 mag. We also required that the source 
have an image m orphology consistent with stellar sources 
(|lrwin et al.ll2007h . Despite our attempts to remove system- 
atic trends in the light curves, we invariably suffered from 
residual correlated red noise, so we modified the detection 
significance statistic, S, from OCCFI T acco rding to the pre- 
scription (equation 4) of I Pont et al.l (|2006h . to obtain S re d- 
Our transit candidates must then have S re d ^ 5 to sur- 
vive, although we note that this is more permissive with 
respect to the limit recommended by Pont and collabora- 
tors (Sred ^ 7). We further discarded any detections with a 
period in the range 0.99 <P< 1.005 days, in order to avoid 
the common ~1 day alias of the ground-based photometric 

surveys. 

Next, as fully described in iBirkbv et al.l (|2012al ). the 
WFCAM ZYJHK single epoch photometry was combined 
with five more optical photometric data points (ugriz 
bands) available for the 19hr field from the Sloan Digi- 
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tizcd Sky Survey archive (SDSS 7 th release, iNashl Il996p 
to create an initi al SEP. The SEDs were fitted with the 
NextGen models jBaraffe et al.|[l99ct ) in order to estimate 
effective temperatures and hence a stellar radius for each 
source. We could then impose a final threshold on the de- 
tected transit depths, by rejecting those that corresponded 
to a planetary radius greater than 2Rj. It is worth noting 
that OCCFIT tends to under-estimate transit depths because 
it does not consider limb-darkening effects and the trape- 
zoidal shape of the transits. Moreover, the NextGen models 
syste matically under-esti mate the temperature of solar-like 
stars |Baraffe et al.l|T998f ). so the initial radius estimates are 
also under-predicted. Hence, the genuine HJ transit events 
are unlikely to be removed in this final threshold cut. As 
a result, our transits detection procedure was conservative. 
Many of the ~3 500 phase-folded light curves which satis- 
fied our criteria were false-positives arising from nights of 
bad data or single bad frames. Others were binary systems, 
folded on half the true orbital period. A more detailed anal- 
ysis of the candidate selection procedure and of advanced 
selecti on steps in the survey can be found in ISipocz et al.l 
l|2012l ). 

The WTS-1 J-band light curve passed all our selection 
criteria and the object (see Table[2]) progressed to the follow- 
ing phases of candidate's confirmation. The descriptions of 
the photometric and spectroscopic data in the next sections 
are organized in order to match the following chronologi- 
cal sequence of analysis. First, optical i'-band photometric 
follow-up of WTS-1 was conduced in order to prove the real 
presence of the transit and check the transit depth consis- 
tency between the two bands. Then, ISIS/WHT intermedi- 
ate resolution spectra enabled us to place strong constrains 
on the velocity amplitude of host star (at the km s^Mevel) 
and therefore to rule out the false-positive eclipsing binaries 
scenarios. Finally we moved to the high-resolution spectro- 
scopic follow-up. The confirmation of the planetary nature 
of WTS-lb came with the RV measurements obtained with 
the HET spectra. 



2.1.3 Broad band photometry 

The WFCAM and SDSS photometric data for WTS-1 are 
reported in Table [3] with other single epoch broad band pho- 
tometric observations. Johnson B V R bands were observed 
for WTS-1 on the night of 6th April 2012 at the Univer- 
sity of Hertfordshire's Bayfordbury Observatory. We used 
a Meade LX200GPS 16- inch f/10 telescope fitted with an 



Table 3. Broad band photometric data of WTS-1 measured 
within the WFCAM (Vega), SDSS (AB), 2MASS (Vega) and 
WISE (Vega) surveys. Johnson magnitudes in the visible are pro- 
vided too (Vega). Effective wavelength ^ e ff (mean wavelength 
weighted by the transmission function of the filter), equivalent 
width (EW) and magnitude arc given for each single pass-band. 
The bands arc sorted by increasing ^ e ff- 



Band 


*e//[A] 


EW[A] 


Mag 


nitude 


SDSS-u 


3546 


558 


18.007 


(±0 


014) 


Johnson-B 


4378 


970 


17.0 


(±0 


1) 


SDSS-g 


4670 


1158 


16.785 


(±0 


004) 


Johnson-V 


5466 


890 


16.5 


(±0 


1) 


SDSS-r 


6156 


1111 


16.434 


(±0 


004) 


Johnson-R 


6696 


2070 


16.1 


(±0 


1) 


SDSS-i 


7471 


1045 


16.249 


(±0 


004) 


UKIDSS-Z 


8817 


879 


15.742 


(±0 


005) 


SDSS-z 


8918 


1124 


16.189 


(±0 


008) 


UKIDSS-Y 


10305 


1007 


15.642 


(±0 


007) 


2MASS-J 


12350 


1624 


15.375 


(±0 


052) 


UKIDSS-J 


12483 


1474 


15.387 


(±0 


005) 


UKIDSS-H 


16313 


2779 


15.103 


(±0 


006) 


2MASS-H 


16620 


2509 


15.187 


(±0 


081) 


2MASS-Ks 


21590 


2619 


15.271 


(±0 


199) 


UKIDSS-K 


22010 


3267 


15.048 


(±0 


009) 


WISE-W1 


34002 


6626 


15.041 


(±0 


044) 


WISE-W2 


46520 


10422 


15.886 


(±0 


157) 



SBIG STL-6303E CCD camera, and integration times of 300 
seconds per band. Images were bias, dark, and flat-field cor- 
rected, and extracted aperture photometry was calibrated 
using three bright reference stars within the image. Photo- 
metric uncertainties combine contributions from the SNR of 
the source (typically ~20) with the scatter in the zero point 
from the calibration stars. The Two Micron All Sky Survey 
f2MASS, ISkrutskie et ahll2006l) and the Wi de-field Infrared 
Survey Explorer (WISE. I Wright et al.ll2010T ) provide further 
NIR data points (J H Ks bands and Wl W2 bands respec- 
tively) . 



2.1.4 INT i'-band data 

In addition to the WFCAM J-band light curve, we observed 
one half transit of the WTS-1 system in the i'-band using 
the W ide Field Camera on the 2.5-m Isaac Newton Tele- 
scope (|McMahon et alj 120011 ) on July 23, 2010. A total of 



© 2012 RAS, MNRAS OOO.ITJITSl 



Inflated HJ in close orbit around a late F-star 5 



Table 4. INT i'-band light curve of WTS-1. The full epoch list, 
which contains 82 entries, will be available electronically. 



HJD 


Normalized 


Error 


-2 400 000 


flux 




55401.3703254 


1.0244 


0.0109 


55401.3809041 


1.0205 


0.0034 


55401.4894461 


0.9936 


0.0022 



82 images, sampling the ingress of the transit, were taken 
with an integration time of 60 seconds. The data were re- 
duced with the CASU INT/WFC data redu ctio n pipeline as 
descri bed in detail bv llrwin fc Lewisl (|200lh and llrwin et al.l 
l|2007h . The pipeline performs a standard CCD reduction, 
including bias correction, trimming of the overscan and non- 
illuminated regions, a non-linearity correction, flat-fielding 
and defringing, followed by astrometric and photometric cal- 
ibration. A master catalogue for the i'-band filter was then 
generated by stacking 20 frames taken under the best con- 
ditions (seeing, sky brightness and transparency) and run- 
ning source detection software on the stacked image. The 
extracted source positions were used to perform variable 
aperture photometry on all of the images, resulting in a 
time-series of differential photometry. 

The final out-of-transit RMS in the WTS-1 i'-band light 
curve is 0.0026 (equivalent to 2.87 mmag) and is used to 
refine the transit model fitting procedure in Section 13.2.11 
The i'-band light curve for WTS-1 is given in Table [4] 

2.2 Spectroscopic data 

2.2.1 ISIS/WHT 

We carried out intermediate-resolution spectroscopy of the 
star WTS-1 over two nights between July 29 - 30, 2010, 
as part of a wider follow-up campaign of the WTS planet 
candidates, using the William Herschel Telescope (WHT) at 
Roque de Los Muchachos, La Palma. We used the single-slit 
Intermediate dispersion Spectrograph and Imaging System 
(ISIS). The red arm with the R1200R grating centred on 
8500 A was employed. We did not use the dichroic during 
the ISIS observations because it can induce systematics and 
up to 10% efficiency losses in the red arm, which we wanted 
to avoid given the relative faintness of our targets. The 
four spectra observed have a wavelength coverage of 8100- 
8900 A. The wavelength range was chosen to be optimal for 
the majority of the targets for our spectroscopic observa- 
tion which were low-mass stars. The slit width was chosen to 
match the approximate seeing at the time of observation giv- 
ing an average spectral resolution R ~ 9000. An additional 
low-resolution spectrum was taken on July 16th, 2010, us- 
ing the ISIS spectrograph with the R158R grating centred on 
6500 A. This spectrum has a resolution (R ~ 1000), a SNR of 
~ 40 and a wider wavelength coverage (5000-9000 A). The 
spectra were processed using the IRAF.CCDPROC0 package 

3 IRAF is distributed by National Astronomy Observatories, 
which is operated by the Association of Universities of Research 
in Astronomy, Inc., under contract to the National Science, USA 



for instrumental signature removal. We optimally extracted 
the spectra and performed wavelength calibration using the 
semi-automatic kpno.doslit package. The dispersion func- 
tion employed in the wavelength calibration was performed 
using CuNe arc lamp spectra taken after each set of expo- 
sures. 



2.2.2 CAFOS/2.2-m Calar Alto 

Two spectra of WTS-1 were obtained with CAFOS at the 
2.2-m telescope at the Calar Alto observatory (as a Direc- 
tors Discretionary Time - DDT) in June, 2011. CAFOS is a 
2kx2k CCD SITe#ld camera at the RC focus, and it was 
equipped with the grism R-100 that gives a dispersion of ~ 
2.0 A/pix and a wavelength coverage from 5850 to 9500 A, 
approximately. Their resolving power is of around R ~ 1900 
at 7500 A, with a SNR ~ 25. The data reduction was per- 
formed following a standard procedure for CCD processing 
and spectra extraction with iraf. The spectra were finally 
averaged in order to increase the SNR. 



2.2.3 KPNO 

A low resolution spectrum of WTS-1 was observed in 
September 2011 with the Ritchey-Chretien Focus Spectro- 
graph at the 4-m telescope at Kitt Peak (Arizona, USA). 
The grism BL-181, which gives a dispersion of ~ 2.8 A/pix, 
was used. Calibration, sky subtraction, wavelength and flux 
calibration were performed following a standard procedure 
for long slit observations using dedicated iraf tasks. The 
ThAr arc lamp and the standard star spectra, employed 
for the wavelength and the flux calibration respectively, 
were taken directly after the science exposure. The mea- 
sured spectrum covers the wavelength range 6000-9000 A 
with SNR~40 at 7500 A and has a resolution of R ~ 1000. 
This is the only flux calibrated spectrum we have available. 



2.2.4 HET 

In the late 2010 and in the second half of the 2011 the star 
WTS-1 was observed during 11 nights with the High Res- 
olution Spectrograph (HRS) housed in the insulated cham- 
ber in the basement of the 9.2-m segmented mir ror Hobby- 
Eberly Telescope (HET, see lRamsev et alj|l998h . The HRS 
is a single channel adaptation of the ESO UVES spectrome- 
ter linked to the corrected prime focus of the HET through 
its fiber- fed instrument as described by iTullI (|l99&f ). It uses 
an R-4 echelle mosaic, which we used with a resolution of 
R=60 000, and a cross-dispersion grating to separate spec- 
tral orders, while the detector is a mosaic of two thinned 
and anti-reflection coated 2K x 4K CCDs. 

One science fiber was used to get the spectrum of the 
target star and two sky fibers were used in order to subtract 
the sky contamination. A couple of ThAr calibration expo- 
sures were taken immediately before and after the science 
exposure. This strategy allowed us to keep under control any 
undesired systematic effect during the observation. Each sci- 
ence observation (except one, see Section 2. 2 \ was split in 
two exposures, of about half an hour each, in order to limit 
the effect of the cosmic rays hits, which can affect the data 
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reduction and analysis steps. Due to the faintness of the 
star, the Iodine gas cell was not used for our observations. 

The data reduction was performed w it h th e 
IRAF.ECHELLE package (|Willmarth fc Barnesl 1 1994 ). 
After the standard calibration of the raw science frame, 
bias-subtraction and flat-fielding, the stellar spectra were 
extracted order by order and the related sky spectrum 
was subtracted. The extracted ThAr spectra were used to 
compute the dispersion functions, which are characterized 
by an RMS of the order of 0.003 A. Consistency between the 
two solutions (computed from the ThAr taken before and 
after the science exposure) was checked for all the nights 
in order to detect possible drifts or any other technical 
hitch that could take place during each run. Successively, 
the spectra were wavelength calibrated using a linear 
interpolation of the two dispersion functions. 

In the subsequent data analysis, custom Matlab pro- 
grams were used. After the continuum estimation and the 
following normalization, the spectra were filtered to remove 
the residual cosmic rays peaks left after the previous filtering 
performed on the raw science frame with the iraf task COS- 
MICRAYS. Comparing the spectra of all the nights for each 
single order, the pixels with higher flux, due to a cosmic hit 
on the detector, were detected and masked. The telluric ab- 
sorption lines present in the redder orders were then removed 
in order to avoid contaminations using a high SNR observed 
spectrum of a white dwarf as template for the telluric lines. 
Finally, the spectra related to the two split science expo- 
sures were averaged to obtain the final set of spectra used in 
the following analysis. A total of 40 orders (18 from the red 
CCD and 22 from the blue one) cover the wavelength range 
4400-6300 A. The spectra have a SNR ~ 8. 



3 ANALYSIS AND RESULTS 
3.1 Stellar parameters 

3.1.1 Spectral energy distribution 

A first characterization of the parent star can be performed 
comparing the shape of the Spectral Energy Distribution 
(SED), constructed from broad band photometric observa- 
tions, with a grid of synthetic theoretical spectra. The data 
relative to the photometric bands, collected in Table|3] were 
analysed with the application VQSA (Virtual Observatory 
SED Analyser, iBavo et alj|200sl . l2012ft . VOSA offers a valu- 
able set of tools for the SED analysis, allowing the estima- 
tion of the stellar parameters. It can be accessed through its 
web-page interface and accepts as input file an ASCII table 
with the following data: source identifier, coordinates of the 
source, distance to the source in parsecs, visual extinction 
(Ay), filter label, observed flux or magnitude and the related 
uncertainty. For our purpose, we tried to estimate only the 
main intrinsic parameters of the star: effective temperature, 
surface gravity and metallicity. The extinction Ay was as- 
sumed as a further free parameter. 

The synthetic photometry is calculated by convolving 
the response curve of the used filter set with the theoreti- 
cal synthetic spectra. Then a statistical test is performed, 
via xt minimization, to estimate which set of synthetic pho- 
tometry best reproduces the observed data. We decided to 
employ the Kurucz ATLAS9 templates set (|Castelli et al.l 




3 5 10 20 30 40 50 60 

Wavelength [1C?A] 



Figure 2. Broad band photometric data of WTS-1: SDSS (filled 
dots), Johnson (empty dots), WFCAM (squares), 2MASS (trian- 
gles) and WISE (stars). The best fitting template (black line) is 
the ATLAS9 Kurucz T cff =6500K, log 3=4.5, [Fe/H]=-0.5 model 
with Ay =0.44. Vertical errorbars correspond to the flux uncer- 
tainties while those along the X-axis represent each band's EW 
(see Table 0. 

11997? ) to fit our photometric data. These templates repro- 
duce the SEDs in the high temperature r egime better than 
the NextGen models (|Baraffe et al.|[l998h . more suitable at 
lower temperatures (< 4500 -ft'). 

In order to speed up the fitting procedure, we restricted 
the range of T cft and log g to 3500-8000 K and 3.0-5.0, re- 
spectively. These constrains in the parameter space did not 
affect the final results as it was checked a posteriori that 
the same results were obtained considering the full avail- 
able range for both parameters. The resulting best fitting 
synthetic template, plotted in Figure [2] with the photo- 
metric data, corresponds to the 7^=6500 K, log <?=4.5 and 
[Fe/H]=-0.5 Kurucz model and Ay=0.44. Uncertainties on 
the parameters were estimated both using xt statistical 
analysis and a bayesian (flat prior) approach. The related 
errors result to be of the order of 250 K, 0.2, 0.5 and 0.07 
for Teff, log g, [Fe/H] and Ay, respectively. 

As it can be seen in Figure^ the WISE W2 data point is 
not consistent with the best fitting model. Firstly, it is worth 
noting that the observed value of W2 — 15.886 (±0.157) is 
below the 5 sigma point source sensitivity expected in the 
VF2-band (> 15.5). The number of single source detections 
used for the VF2-band measurement is also considerably less 
than that of the IVl-band (4 and 19 respectively) increasing 
the uncertainty in the measurement. Finally, the poor an- 
gular resolution of WISE in the W2-band (6.4") could add 
further imprecisions to the final measured flux, especially in 
a field as crowded as the WTS 19hrs field. For these reasons, 
the WISE W2 data point was not considered in the fitting 
procedure. 

Once the magnitude values were corrected for the in- 
terstellar absorption according to the best fitting value 
(Ay=0.44±0.07), colour - temperature relations were used 
to further check the effective temperature and spectral type 
of the host star. From the SDSS g and r magnitudes, we 
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Figure 3. Upper row: comparison of a degraded spectrum of the Sun (black) with a synthetic spectrum (red) computed with T o ff=5750K, 
log 3=4.5 and [Fe/H]=0.0. Lower row: comparison of the WTS-1 spectrum (black) with different synthetic spectra (colours). The 
comparison between the observed WTS-1 spectrum and the synthetic models took into account the differences of the core of the lines 
shown in the upper row plots. From this analysis, the best fitting model is the one with T c ff=6250K, log g=4.5 and [Fe/H]=0.0. 



obtained a value of (B-V) =-0.43±0.04 (| Jester et all 120051 ) 
which imply r e fT=6300±600 K assumi ng log <?=4.4 and 
[Fe/Hl=-0.5 JSekiguchi fc Fukugita!l2000ri . Following the ap- 
pendix B of lCollier Cameron et alT(|2007l ), the 2MASS (J - 
H) index of 0.23 ±0.09 leads t o a va lue of T e fr=6200±400 K 
while Table 3 of lCovev et all l|2007t) suggests the host star 
to be a late-F considering different colour indices at once. 
These results are all compatible within the uncertainties. 



3.1.2 Spectroscopic analysis 

The spectroscopic spectral type determination was done 
firstly by comparing the spectrum observed with CAFOS 
with a set of spectra of template stars. Stars of different 
spectral types, uniformly spanning the F5 to G2 range, were 
observed with the same instrumental setting. Since the ob- 
served spectrum has a relatively low SNR, we focused the 
analysis on the strongest features present which are the 
Ha (6562.8 A) and the Call triplet (8498.02 A, 8542.09 A, 
8662.14 A). The best match was obtained, via minimization 
of the RMS of the difference between the WTS-1 and tem- 
plate star spectra, with the spectrum of an F6V star with 
solar metallicity. 

Afterwards, we tried to estimate the stellar parame- 
ters comparing the observed spectrum with a simulated 
spectrum with known parameters. For that aim, we used 
a library of high resolution synthetic stellar spe ctra by 
ICoelho et all l|2005h . created by the PFANT code i|Barbuvl 



1 19821 ; ICavrel et al]|l99ll ; iBarbuv et a l1 l2003 j) that computes 
a synthetic spectrum assuming local thermodynamic equilib- 
rium (LTE). The synthetic spect ra were achieved us i ng th e 
model atmospheres presented bv lCastelli fc Kurucj l|2003h . 
Since the core of these lines are strongly affected by cu- 
mulative effects of the chromosphere, non-LTE (local ther- 
modynamical equilibrium) and inhomogeneity of velocity 
fields, we firstly compared a spectrum of the S un observed 
with HIRES spectrograph at the Keck telescope (|Vogt et alj 
1 1994 ). The spectrum of the Sun was degraded to lower reso- 
lution and resampled to match our CAFOS spectrum speci- 
fications to see how such effects appear at this resolution and 
how different the solar spectrum is from a synthetic spec- 
trum from the library by Coelho and collaborators. Look- 
ing at the upper plots in Figure [3] we concluded that the 
cores of the lines in the simulated spectrum are systemati- 
cally higher than those in the observed spectrum of the Sun. 
Nevertheless, the Call triplet line at 8498.02 A seems to be 
less affected by the above-mentioned problems. Consider- 
ing these differences between central depth of the observed 
and simulated spectra in the best fitting procedure, we esti- 
mated that the synthetic model with T o ff=6250 K, log <7=4.5 
and [Fe/H]=0.0 best reproduces the observed spectrum of 
WTS-1. The expected uncertainties are of the order of the 
step size of the used library (<5T , o fr=250 K, <51og g=0.5 and 
5[Fe/H]=0.5). 

The high resolution HET spectra were employed to at- 
tempt a more detailed spectroscopic analysis of the host 
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ISIS spectrum 
Teff=6000; logg=4.4 
Teff=6250; logg=4.4 
Teff=6500; logg=4.4 
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Figure 4. Intermediate resolution ISIS spectrum (black line) and 
high resolution HET spectrum (grey line) of the Ha line of WTS- 
1 in comparison with synthetic spectra calculated with different 
effective temperatures of 6000 K (blue dash-dotted line), 6250 K 
(dashed line; finally adopted T c g), and 6500 K (blue dotted line). 
The synthetic spectra were convolved in order to match the res- 
olution of the ISIS spectrum. 



star. The spectra observed each single night were stacked 
together obtaining a final spectrum with a SNR of about 
12, calculated over a 1 A region at 5000 A. To compute 
model atmospheres, LLm odels stellar model atmosphere 
code (|Shulvak et all |2004| ) was used. For all the calcula- 
tions, LTE and plane-pa rallel geometry were assumed. We 
used t he VALD database ( Piskunov et al Jl995l :l iKupka et al.l 
1 19991 ; IRvabchikova et alj ll999f T as a source of atomic line 
parameters for opacity calculations with the LLmodels 
code. Finally, convec t ion w as implemented according to the 
ICanuto fc Mazzitellil |l99ll ) model of convection. 

The parameter determination and abundance analysis 
were performed iteratively, self-consistently recalculating a 
new model atmosphere any time one of the parameters, in- 
cluding the abundances, changed. As a starting point, we 
adopted the parameters derived from the CAFOS spectrum. 
We performed the atmospheric parameter determination by 
imposing the iron excitation and ionization equilibria mak- 
ing use of equivalent widths measured for all available un- 
blended and weakly blended lines. We converted the equiva- 
lent width of each line into a n abundance value with a mod- 
ified version (|Tsvmballll99r3 ) of the WIDTH9 code (jKurucj 
1 19931 ). Unfortunately, the faintness of the observed star, cou- 
pled with the calibration process (including the sky subtrac- 
tion), led to a distortion of the stronger lines, weakening 
their cores. For this reason, our analysis took into account 
only the measurable weak lines, making therefore impossi- 
ble a determination of the microturbulenc e velocity (vmic), 
which we fixed at a value of 0.85 km s _1 ()Valenti fc Fishen 
120051 ). With the fixed v m ic, we imposed the Fe excitation and 
ionization equilibria, which led us to an effective tempera- 
ture of 6000±400K and a surface gravity of 4.3±0.4. Im- 
posing the ionization equilibrium we took into account the 
expec ted non- LTE effects for Fe I (~0.05 dex, iMashonkinal 
l201ll ). while for Fe 11 non-LTE effects are negligible. 

In this temperature regime, the ionization equilibrium 



is sensitive to both T e s and log g variations, therefore it is 
important to simultaneously and independently further con- 
strain temperature and/or gravity. For this reason we looked 
at the Ha line to further constrain T c g, as in this temper- 
ature r egime Ha is sensitive p rimarily to temperature vari- 
ations (|Fuhrmann et al]|l993h . We did t his by fitting syn - 
thetic spectra, calculated with SYNTH3 |Kochukhovll2007r i . 
to the Ha line profile observed in the HET high resolution 
and in the ISIS/WHT low resolution spectra. As the Ha line 
of the HET spectrum was also affected by the before men- 
tioned reduction problems, only the wings of the line were 
considered. Although we could calculate synthetic line pro- 
files of Ha on the basis of atmospheric models with any T e g , 
because of the low SNR of our spectra, we performed the line 
profile fit ting on the basis of mod els with a temperature step 
of 100 K (jFuhrmann et al.lll993l . small variations in gravity 
and metallicity are negligible). From the fit of the Ha line we 
obtained a best fitting T c ff of 6100±400 K. Furt her details on 
method, codes and techniques ca n be found in Fossati et al.l 
|2009f ). IRvabchikova et all (|2009l ). iFossati et al.1 |201lf ) and 



references therein. Figure [4] shows the Ha line profile ob- 
served with ISIS and HET in comparison with synthetic 
spectra calculated assuming a reduced set of stellar parame- 
ters. On the basis of the previous analysis, we finally adopted 
Teff=6250±250 K, log g=4.4±0.1. With this set of param- 
eters and the equivalent widths measured with the HET 
spectrum, a final metallicity of -0.5±0.5 dex dex was derived, 
where the uncertainty takes into account the internal scatter 
and the uncertainty on the atmospheric parameters. By fit- 
ting synthetic spectra, calculated with the final atmospheric 
parameters and abundances, to the observation, we derived 
a projected rotational velocity v sm(i)=7±2km s~ . 

The HET spectrum allowed us to measured the atmo- 
spheric lithium abundance from the Li I line at ~6707A. 
Lithium abundance is important as it can constrain the 
age of the star (see Section I3,1.3[) . As this line presents 
a strong hyperfine structure and is slightly blended by a 
nearby Fel line, we measured the Li I abundance by means 
of spectral synthesis, instead of equivalent widths. By adopt- 
in g the meteoritic /terrestr ial isotopic ratio Li6/Li7=0.08 
bv lRosman fc Tavlorl (|l998T ). we derived log N(Li)=2.5±0.4 
(corresponding to an equivalent width of 41.12±24.40 mA), 
where the uncertainty takes into account the uncertainty on 
the atmospheric parameters, T e ff in particular (see Figure[5]). 

The low resolution spectrum obtained at the KPNO 
observatory, being the only flux calibrated spectrum, was 
compared to the Kurucz ATLAS9 templates set (the same 
employed in the SED analysis in Section 13.1.1(1 . The lim- 
ited wavelength range covered by the observed spectrum 
(3000 A) , allowed to achieve usable results fitting only one 
parameter of the template spectra. We therefore decided to 
leave T e g as a free parameter, fixing the other quantities 
to log 5=4.5, [Fe/H]=-0.5 and Av=0.44. We concluded that 
the temperature range 6000-6500 K brackets the T e ft with 
95% confidence level. 



3.1.3 Properties of the host star WTS-1 

The atmospheric parameters of the star WTS-1 were com- 
puted combining the results coming from the analysis de- 
scribed in the previous sections. We finally obtained an effec- 
tive temperature of 6250±200 K, a surface gravity of 4.4±0.1 
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Table 5. Properties of the WTS-1 host star. 



Parameter 



Value 



6707.5 6708.0 6708.5 

wavelength (A) 

Figure 5. High-resolution HET stacked spectrum (black full line) 
of the WTS-1 Li I line at ~6707A in comparison with three syn- 
thetic spectra calculated with the final adopted stellar parameters 
and lithium abundances (log N(Li)) of 2.5 (red dashed line), 2.9 
(blue dotted line), and 2.1 (blue dash-dotted line). The vertical 
line indicates the position of the centre of the multiplet. 



and a metallicity range —0.5 -f- 0.0 dex. As described in Sec- 
tion ^, 1.2l the faintness of the star and the reduction process 
led to a distortion of the stronger lines in the HET high res- 
olution spectrum, reducing the number of reliable lines em- 
ployed in the measure of the metallicity. For these reasons, 
the relative uncertainty on the metallicity is larger with re- 
spect to those of the other parameters. Further observations 
would be needed to pin down the exact value of the star 
metallicity. 

In order to determine the parameters of the stellar com- 
panion, mass and radius of the host star must be known. 
The fit of the transit in the light curves provides important 
constrains on the mean stellar density (see Section [3.2. 1 [) . 
Joining this quantity to the effective temperature, we could 
place WTS-1 in the modified p a vs. T a g H-R diagram and 
compar e its position with ev olutionary tracks and isochrones 
models (|Girardi et al.ll2000l ) in order to estimate stellar mass 
and age. In this way, we estimated the stellar mass to be 
1.2±0.1 Mq and the age of the system to range between 
200 Myr and 4.5 Gyr. 

Further constrains on the stellar age can be fixed consid- 
ering the measured Li I line abundance. Depletion of lithium 
in stars hotter than the Sun is thought to be due to a not 
yet clearly identified slow mixing process during the main- 
sequence evolution, because those stars d o not experience 
pre-main sequence depletion (|Martfnlll997l ). Comparison of 
the lithium abundance of WTS-1 (log N(Li)=2.5±0.4) with 
those of open clusters rise the lower limit on the age to 
600 Myr, because younger cluster s do not show lithium de - 
pletion in their late-F members jSestito fc Randichll2005h . 



On the other hand, it is not feasible to derive an age con- 
strain from the WTS-1 rotation rate (v sm(i)=7±2km s -1 ) 
because stars with sp ectral types earlie r than F8 show no 
age-rotation relation (| Wolff et al.lll986l ) and thus they are 
left out fro m the formula tion of the spin down rate of low- 
mass stars (|Stepienlll98sh . 

The true space motion knowledge of WTS-1 allows to 
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determine to which component of our Galaxy it belongs. 
In order to compute the U, V and W components of the 
motion, the following set of quantities was required: dis- 
tance, systemic velocity (from the RV fit, see Section 3.2.21) . 
prop er motion (from SDSS 7 th release, iMunn et al 2004, 
l2008h and coordinates of the star. The distance to the ob- 
served system was estimated according to the extinction, 
fitted in the SED analysis (Ay =0. 44 ± 0.07), and a mode l 
of dust distribution in the galaxy (jAmores fc Lepinej I2005L 
axis-symmetric model). The UVW values and their errors 
are ca lculated using the method in Ijohnson fc Soderbloml 
(jl987h . with respect to the Sun (heliocentric) and in a left- 
handed coordinate system, so that they are positive away 
from the Galactic centre, Galactic rotation and the North 
Galactic Pole respectively. All the quantities here discussed 
are listed in Table [5] Considering the uncertainties on the 
derived quantities, maily affected by the error on the dis- 
tance, the host star is consistent with both the definitions of 
Galactic young-old disk and young disk populations (metal- 
licity between -0.5 and Odex, solar- metallicity respectively, 
lLeggettJ[l99l ). But we could assess that WTS-1 is not a halo 
member. 



Combining p a and M 3 , a value of 1.15to'i2 ^0 was com- 
puted for the stellar radius. The same result was obtained 
considering the stellar mass and the surface gravity mea- 
sured from the spectroscopic analysis. Scaling by the dis- 
tance the apparent V magnitude calcula ted from the SDSS 
g and r magnitudes (jjester et al.l 120051 ) , we computed an 
absolute V magnitude of 3.55lQ'gg. This value and all the 
other derived stellar parameters are consistent with each 
other and with the typical quantities e xpected for an F6-8 
main-sequence star (I Torres et al.|[2009l ). They are collected 
in Table H 
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Table 6. Quadratic limb-darkening coefficients used for the tran- 
sit fitting, for a star with effective temperature T e ff=6250K, 
surface gravity log .9= 4.4 and mctallicity [Fe/H]=-0.5dex 
llClaret fe Bloemedl20nh . 



filter 


71 


72 


J 


0.14148 


0.24832 


i' 


0.25674 


0.26298 



3.2 Planetary parameters 

3.2.1 Transit fit 

The light curves in J- a nd i'-band were fitted with ana- 
lytic models presented by llvlandel fc Ago! (|2002u . We used 
quadratic limb-darkening coefficients for a star with effec- 
tive temperature T e fr=6250K, surface gravity log <?=4.4 and 
metallicity [Fe/H]=-0.5 dex, calculated as linear interpola- 
tions in T e g, log q and [Fe/H] of the values tabulated in 
IClaret fc Bloemenl (|201ll ). We use their table derived from 
ATLAS atmospheric models using the flux conservation 
method (FCM) which gave a slightly better fit than the ones 
derived using the least-squares method (LSM). The values 
of the limb-darkening coefficients we used in our fitting are 
given in Table [6] Scaling factors were applied to the error 
values of the J- and i'-band light curves (0.94 and 0.9 re- 
spectively) in order to achieve a reduced \ 2 °f the constant 
out-of-transit part equal to 1. 

Using a simultaneous fit to both light curves, we fitted 
the period P, the time of the central transit to, the radius 
ratio R P /R B , the mean stellar density, p s = M B /R B in solar 
units and the impact parameter /3i mpa ct in units of R B . The 
light curves and the model fit are shown in Figures [6] and [7] 
while the resulting parameters are listed in Table [7] 

The errors were calculated using a multi-dimensional 
grid on which we search for extreme grid points with Ax 2 =l 
when varying one parameter and simultaneously minimizing 
over the others. Figure [8] shows the correlations between the 
parameters of the WTS-1 system derived from the simulta- 
neous fit of the J- and i'-band light curves. Considering the 
solar metallicity scenario, with different limb-darkening co- 
efficients, the change in the final fitting parameters is smaller 
than 1% of their uncertainties. Note that the combined fit as- 
sumes a fixed radius ratio although in a hydrogen-rich atmo- 
sphere, molecular absorption and scattering processes could 
result in different radius ratios in each band (an attempt 
to detect such varia tions has recently been undertaken by 
Ide Mooii et~al1l2012u . In our case, the uncertainties are too 
large to see this effect in the light curves and the assump- 
tion of a fixed radius ratio is a good approximation. The 
estimated stellar density of the host star (0.791q p SU n) is 
consistent with the expected val ue based on its spectral type 
(|Seager fc Mallen-Ornelasll2003l ). The noise in the data did 
not allow a secondary transit detection. 

Subsequently, we searched for further periodic signals in 
the light curve after the removal of the data points related to 
the transit. No significant signals were detected in the Lomb- 
Scargle periodogram up to a period of 400 days. Since WTS- 
1 is a late F-star, there are not many spots on the surface 
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Figure 6. WFCAM J-band light curve data of WTS-1. Upper 
panel: whole set of the folded WFCAM J-band data points. 
Middle panel: folded photometric data centred in the transit and 
best model (red line) fitted in combination with the INT i'-band 
data. Lower panel: residuals of the best fit. 
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Figure 7. INT i'-band light curve data of WTS-1. Upper panel: 
photometric data and best model (red line) fitted in combination 
with the WFCAM J-band data. Lower panel: residuals of the 
best fit. 



and they do not live long enough to produce a stable signal 
over a timescale of several years. 



3.2.2 Radial velocity 

One of the most pernicious transit mimics in the WTS are 
eclipsing binaries. On one hand, a transit can be mimicked 
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Figure 8. Correlation plots of the quantities derived from the simultaneous fit of the transit in the J and i'-band light curves. The x 2 
minima are indicated by crosses while the different tones of grey correspond to the 68%, 95% and 99% confidence level (darker to lighter 
respectively). The other couples of parameters do not show significant correlations. 



Table 7. Fitted parameters of the WTS-f system as determined 
from the simultaneous fit of the J- and i'-band light curves. Scal- 
ing factors to the uncertainties of the J and %' data points (0.94 
and 0.9 respectively) were applied in order to achieve a reduced 
X 2 =l in the constant out-of-transit part of the light curves. 



Parameter 



Value 



Porb 

to 

Rj,/R s 

Ps 

Pimpact 



3.352059 
2 454 318.7^ 
0.1328 



+ 1.2x10 
1.4x10 

0.0043 
0.0036 



days 
HJD 



j+0.0032 
-0.0035 



7Q+ 31 „ 
' y _n is Psun 



0.69 



-0.18 
+0.05 
0.09 



by an eclipsing binary that is blended with foreground or 
background star. On the other hand, grazing eclipsing bi- 
naries with near-equal radius stars also have shallow, near- 
equal depth eclipses that can phase-fold into transit-like sig- 
nals at half the binary orbital period. In order to rule out the 
eclipsing binaries scenarios, the RV variation of the WTS-1 
system were first measured using the ISIS/WHT intermedi- 
ate resolution spectra. Eclipsing binaries systems typically 
show RV amplitudes of tens of km s _1 , while the measured 
RVs were all consistent with a flat trend within the RV un- 
certainties of ~ 1 km s _1 . 

Afterwards, we analysed the high-resolution HET spec- 
tra in order to accurately investigate the properties of the 
sub-stellar companion of WTS-1. The spectra related to each 
single night were cross-correlated using the iraf.rv.fxcorr 
task with the synthetic spectrum of a star with T e ff =6250 K, 
log g~4:A and [Fe/H]=-0.5. Changes in the effective temper- 
ature of the synthetic templates, even of the order of sev- 
eral hundreds of Kelvin, cause variations of the measured 
RV values smaller than the statistical uncertainties due to 
noise in the spectra (suggesting the absence of contami- 
nation of back/foreground stars of different spectral type). 
Even smaller variations occur changing surface gravity and 
metallicity of the template. The 40 single RV values, each 
one coming from a different order, were used to compute 
the RVs and related uncertainty at each epoch of the obser- 
vations. Resampling statistical tools were used in order to 




0.4 0.6 
Phase 



Figure 9. Top panel: RVs values measured with the high- 
resolution HET spectra of WTS-1 as a function of the orbital 
phase. Black dots and blue triangles refers to observations per- 
formed in 2010 and 2011 respectively. The data point at <p = 0.33, 
empty dot, was excluded from the fitting procedure (see text for 
details). Best-fit circular orbit model (x^=1.45) and la uncer- 
tainty of the semi- amplitude {K *=479±34 m s" 1 ) are indicated 
with solid and dashed red lines respectively. The black dotted line 
refers to the fitted radial systemic velocity (7=-l 714±35m s _1 ). 
Zero phase corresponds to the mid-transit time. Middle panel: 
Phase folded O-C residuals from the best-fit. The residual scatter 
is of the order of ~ 150 m s _1 , consistent with the RVs uncer- 
tainties. Bottom panel: Bisector spans. No significant deviations 
from zero and no correlation with the RVs were found. 
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Table 8. Radial velocities and bisector spans measurements for 
WTS-1 obtained by HET spectra. The phases were computed 
from the epochs of the observations expressed in Julian date and 
using the P and to values found with the transit fit. 



Table 9. Properties of the new extrasolar planet WTS-lb. 



HJD 


Phase 


RV 




BS 


-2 400 000 




[km 


s 


- 1 ] 


[km s" 1 ] 


55477.676 


0.74 


-1.46 


± 


0.12 


-0.31 ±0.40 


55479.666 


0.33 


-1.80 


± 


0.10 


-0.77 ±0.39 


55499.608 


0.28 


-2.03 


± 


0.10 


0.76 ±0.51 


55513.586 


0.45 


-2.06 


± 


0.17 


0.87 ±0.63 


55522.556 


0.13 


-2.09 


± 


0.44 


-0.39 ±0.35 


55523.542 


0.42 


-1.70 


± 


0.16 


-0.13 ±0.44 


55742.729 


0.26 


-1.97 


± 


0.15 


0.27 ±0.68 


55760.673 


0.79 


-1.11 


± 


0.12 


-0.47 ±0.59 


55782.837 


0.22 


-2.23 


± 


0.07 


0.41 ±0.39 


55824.722 


0.25 


-2.31 


± 


0.06 


0.16 ±0.56 


55849.650 


0.75 


-1.19 


± 


0.06 


0.18 ±0.31 



Parameter 



M p 
Hp 

Prot 



inc 

impact 

to 
Pp 

Ti a 
eq 



Value 



4.01±0.35 Mj 

i-48±8:iS r j 

3.352057+ 1 1 ;^J°: 5 6 
0.047±0.001 AU 
< 0.1 (C.L.= 95%) 
85. 5" 1 



69+ 05 

2 454 318.7472+°'°°3^ HJD 
1.61±0.56 g cm,- 3 , 1.21±0.42 pj 
1500±100K 



a We assumed a Bold albedo Ag=0 and re-irradiating 
fraction F= 1; 



better estimate mean value, standard deviation and possi- 
ble bias in the sample of measured RVs. Finally, the mea- 
sured RVs were corrected for the Earth orbital movements 
and reduced to the heliocentric rest-of-frame. The phase val- 
ues (f> were computed from the epochs of the observations, 
expressed in Julian date, and using the extremely well de- 
termined P and to values (relative uncertainties are of the 
order of 1CF 6 and 10 -9 respectively) obtained from the pho- 
tometric fit (see Table [7}. 

The data, listed in Table [8] were then fitted with a 
simple two parameters sinusoid of the form: 

RV = 7 + K+ain{2Tt4>) (1) 

where is the RV semi-amplitude of the host star 
and 7 is the systemic velocity of the system. Thanks to the 
acquisition of two ThAr calibration exposures (before and 
after the science exposure, see Section I2.2.4|l , we detected 
a small drift between the ThAr lines occurring during the 
science exposures on November 22, 2010. In the presence of a 
suspected systematic trend, which could affect the measured 
RV value, we performed the fitting procedure excluding the 
data point related to that night (at (j> = 0.33). The larger RV 
error of the data point at <j> = 0.13 is due to the integration 
time (half hour) of the science frame which is shorter than 
those of all the other data points (one hour) . The best fitting 
model (%^=1.45) was obtained for _R"*=479±34 m s _1 and 
7=-1714±35m s _1 and is plotted in Figure |9] with the RV 
data. We imposed the orbit to be circular as the eccentricity 
was compatibl e with zero when a K eplerian orbit fit was 
performed (see lAnderson et al.l (|2012l ) for a discussion of the 
rationale for this). In accordance to the RV uncertainties, a 
relatively loose upper limit can be plaved on the eccentricity 
(e < 0.1, C.L - 95%). The fitted RV semi-amplitude implies 
a planet mass of 4.01±0.35 Mj assuming a host star mass 
of 1.2±0.1 Mq (see Section [3.1,3p . The uncertainty on the 
planet mass in mainly driven by the uncertainty on the mass 
of the host star. As can be seen in Figure [9] the RVs related 
to observations performed in late 2010 and in the second 
half of 2011 are consistent, showing no significant long term 
trends in our measurements. 



The HET spectra were employed also to investigate 
the possibility that the measured RVs are not due to true 
Dopplcr motion in response to the presence of a planetary 
companion. Similar RV variations can rise in case of dis- 
tortions i n the line profiles d ue to stellar atmosphere os- 
cillations ([Queloz et al.l l200l|) . To assert that this is not 
our case, we used the same cross-correlation profiles pro- 
duced previously for the RV calculation to compute the bi- 
sector spa ns (BS herea f ter) w hich are reported in Table [5] 
Following [Torres et alj (|2005h . we measured the difference 
between the bisector values at the top and at the bot- 
tom of the correlation function for the different observation 
epochs. In case of contaminations, we would have expected 
to measure BS values consistently different f rom zero and 
a strong correlation with the measured RVs (| Queloz et al.l 
l200ll ; iMandushev et al.lliooH '). As it can be seen in the bot- 
tom panel of Figure [9] the measured BS do not show signif- 
icant deviation from zero within the uncertainties. No cor- 
relation was detected between the BS and the RV values. In 
this way, contaminations that could mimic the effect of the 
presence of a planet were ruled out. 



4 DISCUSSION AND CONCLUSIONS 

In this paper we announce the discovery of a new transit- 
ing extrasolar planet, WTS-lb, the first detected by the 
UKIRT / WFC AM Transit Survey. The parameters of the 
planet are collected in Table [9] WTS-lb is a ~4 Mj planet 
orbiting in 3.35 days a late F-star with possibly slightly sub- 
solar metallicity. With a radius of 1.49+°' Jg Rj, it is lo- 
cated in the upper part of the mass - radius diagram of 
the known extrasolar planets in the mass range 3-5 Mj 
(see Figure I10[) . The parameters of the other planets are 
taken from www. exoplanet. eu at the time of the publication 
of this work. Planets with only an upper limit on the mass 
and/or on the radius are not shown. It is worth noting that 
only a cut-off of the transit depth, different for each survey, 
could act as a selection effect against the detection of plan- 
ets in this upper portion of the diagram. Larger planetary 
radii imply a deeper transit feature in the light curves and 
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Figure 10. Mass - Radius diagram of the known planets with 
a mass in the range 3-5 Mj (black dots). Labels with the related 
planet name are shown for an easier identification. Planets with 
only an upper limit on mass and/or radius are not shown. The 
blue dashed lines represent the iso-density curves. The green dot- 
dashed lines indicate the plan etary radii at different ages accord- 
ingly to lFortnev et all l|2007l ) (see text for details). Masses, radii 
of the planets are taken from www.exoplanet.eu at the time of 
the publication of this work, while the related uncertainties were 
found in the refereed publication. WTS-lb is shown in red. 



thus, within this mass range, larger object are more easily 
detectable. The properties of WTS-lb, as well as those of 
the other two pla nets present in the upper part of the dia- 
;ram, CoRoT-2b dAlonso et al.ll200ct ) and OGLE2-TR-L9b 
Snellen et alll2009f h~ are not explained within standard for- 
mation and e volution models of isolated gas giant planets 
(|Guillodl2005l ). 

The radius anomaly is at the ~ 2a level considering the 
stellar irradiation that retards the contraction of the planets, 
the distance o f the planet from t he host star and the age 
of the planet (jFortnev et al.l 120071). The models of Fortney 
and collaborators predict indeed a radius of 1.2 Rj for a 
600 Myr-old planet (see their Figure 5, 3 Mj and 0.045 AU 
model) . This radius estimate is an upper limit as 600 Myr is 
the lower limit on the age of the WTS-1 system due to the Li 
abundance (see Section r3.1.3[) . The radius trend shown in the 
figure would suggest an age for WTS-lb less than 10 Myr. 
The same significance on the radius anomaly is obtained 
considering empirical relationships coming from the fit of the 
observed radii as a function of the physical properties of the 
star-planet syste m, such as mass, eq uilibrium temperature 
and tidal heating (|Enoch et alj|20l3 . eq. 10). In any case, a 
rapid migration of WTS-lb inward to the highly- irradiated 
domain after its formation seems required. 

Surface day/night temperature gradients due to the 
strong incident irradiation, are likely to generate strong 
wind activity throug h the planet atmosphere. Recently, 
IWu fc Lithwickl (|2012T) showed how the Ohmic heat ing pro- 
posal ( Batygin fc Stevensonll201ol ; IPerna et al.ll201oT ) can ef- 
fectively bring energy in the interior of the planet and slow 



down the cooling contraction of a HJ even on timescales 
of several Gyr: a surface wind blowing across the planetary 
magnetic field acts as a ba ttery that rises Ohmic dis sipa- 
tion in the deeper layers. In lHuang fc Cummind |2012l ). the 
Ohmic dissipation in HJs is treated decoupling the interior 
of the planet and the wind zone. In this scenario, the radius 
evolution for an irradiated HJ planet (see their Figure 9, 
3 Mj) leads to a value consistent with our observation up to 
3 Gyr. 

Accordingly to iFortnev et all |2008l ), the incident flux 
(the amount of energy from the host star irradiation, per 
unit of time and surface, that heats the surface of the planet) 
computed for WTS-lb (1.12±0.26 ■ 10 9 erg s" 1 cm" 2 ) as- 
signs it to the so called pM class of HJs. This classification 
considers the day-side atmospheres of the highly-irradiated 
HJs that are somewhat analogous to the M- and L-type 
dwarfs. In particular, the predictions of equilibrium chem- 
istry for pM planet atmospheres are similar to M-dwarf 
stars, where absorpt ion by TiO, VO, H2O, and CO is promi- 
nent (jLodd crs 2002). Planets in this class are warmer than 
required for condensation of titanium (Ti)- and vanadium 
(V)-bearing compounds and will possess a temperature in- 
version (which could le ad to a smal ler inflation due to Ohmic 
heating accordingly to lHend (|2012l )) d ue to absorptio n of in - 
cident flux by TiO and VO molecules. IFortnev et"ai] (|2008l ) 
propose that these planets will have large day /night effective 
temperature contrasts and an anomalous brightness in sec- 
ondary eclipse at mid-infrared wavelengths. Unfortunately, 
the SNR in the J-band light curve, due to the faintness of 
the parent star WTS-1, is not high enough for such kind of 
detection (see Section \3. 2. 1[) . 

To conclude, the discovery of WTS-lb demonstrates the 
capability of WTS to find planets, even if it operates in a 
back-up mode during dead time on a queue-schedule tele- 
scope and despite of the somewhat randomised observing 
strategy. Moreover, WTS-lb is an inflated HJ orbiting a 
late F-star even if the project is de signed to search for ex- 
trasolar planets hosted by M-dwarfs. lBirkbv et al.|[2012bl will 
present the second WTS detection, WTS-2b, a Jupiter-like 
planet around a cool K-star. 
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